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1  Introduction 

At  present,  different  types  of  pulsed  high  power  microwave  (HPM)  oscillators  and 
amplifiers  are  being  researched  in  laboratories,  as  well  as  being  manufactured  by  industry. 
They  are  used  in  communication  systems,  radars,  electromagnetic  testing  facilities,  scientific 
research,  and  military  projects.  In  the  experiments  with  relativistic  electron  beams,  a  circular 
waveguide  with  a  horn  antenna  is  often  used  providing  radiation  of  generated  pulse  from  the 
HPM  facility  into  free  space.  Systems  measuring  generated  power  are  usually  located  at  some 
distance  from  the  radiating  antenna.  They  consist  of  the  receiving  antenna  connected  to  the 
device  enabling  the  measurement  of  microwave  pulses. 

To-date  the  main  device  being  used  for  the  measurement  of  the  HPM  pulse  power  is  a 
semiconductor  diode.  It  is  able  to  measure  microwave  pulse  power  up  to  a  few  hundreds  of 
mW  [1].  Therefore,  when  it  is  used  to  measure  HPM  pulses,  the  initial  pulse  has  to  be 
strongly  attenuated.  On  the  one  hand,  large  attenuation  of  the  microwave  pulse  results  in  the 
decrease  of  the  measurement  accuracy.  In  addition,  the  size  and  weight  of  the  measurement 
system  increases,  and  its  control  becomes  complicated.  On  the  other  hand,  when  the  measured 
pulsed  microwave  power  is  of  the  order  of  mW,  the  diode  outputs  a  DC  pulse  with  amplitude 
of  the  order  of  mV.  To  measure  so  small  DC  signal  in  the  presence  of  the  stray  pick-up  and 
electromagnetic  interference  that  are  typical  in  the  environment  of  HPM  sources,  a  shielded 
chamber  with  a  measurement  equipment  situated  at  a  respectful  distance  from  the  HPM 
source  is  usually  employed.  The  installation  of  the  shielded  chamber  confines  the  flexibility 
of  the  measurement  system  and  increases  its  price. 

In  Microwave  Laboratory  of  Semiconductor  Physics  Institute,  we  have  developed  an 
alternative  device  for  the  HPM  pulse  measurement  -  a  resistive  sensor  (RS),  performance  of 
which  is  based  on  the  electron  heating  effect  in  semiconductors.  A  sensing  element  (SE)  of 
the  RS  is  inserted  into  transmission  line  -  usually  the  rectangular  waveguide.  The  electric 
field  of  the  microwave  pulse  traveling  in  the  waveguide  heats  electrons  in  the  SE  its  resistance 
increases  and  by  measuring  this  resistance  change  the  microwave  pulse  power  in  the 
transmission  line  is  determined.  The  main  advantages  of  the  RS  over  the  diode  are  the 
following:  it  can  measure  HPM  pulses  directly,  produces  high  output  signal,  is  overload 
resistant,  and  demonstrates  very  good  long  term  stability  [2].  The  resistive  sensors  developed 
in  our  laboratory  for  rectangular  waveguides  are  used  in  laboratories  dealing  with  HPM  pulses 
worldwide.  Unfortunately,  the  RS  for  the  circular  waveguide  has  not  been  designed  yet. 

Having  in  mind  the  drawback  of  the  measurement  of  HPM  pulses  in  facilities 
employing  circular  waveguides  in  the  output,  our  activities  under  the  current  project  will  be 
concentrated  on  the  elucidation  of  peculiarities  of  the  interaction  of  the  electromagnetic  wave 
propagating  in  the  circular  waveguide  with  a  semiconductor  obstacle  placed  on  the  waveguide 
wall.  Our  task  is  to  find  such  parameters  of  the  obstacle  that  it  can  serve  as  the  prototype  of 
the  SE  for  the  circular  waveguide  RS.  Considering  the  semiconductor  obstacle  from  this  point 
of  view,  the  main  requirements  for  it  can  be  formulated.  First,  the  SE  should  not  insert 
considerable  reflection  in  the  waveguide,  so  the  value  of  VSWR  has  been  set  at  <  1.2.  Second, 
the  RS  should  be  able  to  measure  nanosecond-duration  HPM  pulses;  therefore,  the  DC 
resistance  of  the  RS  should  not  exceed  50  Q.  Third,  the  shape  of  the  SE  should  be  taken  as  a 
plate  that  is  important  for  heat  transfer  from  the  SE.  Finally,  the  flat  frequency  response  of  the 
RS  in  the  waveguide’ s  frequency  band  is  preferable. 

In  the  interim  report  [3],  the  investigations  of  the  interaction  of  the  semiconductor 
obstacle  with  (TE0i)  H0i  mode  electromagnetic  wave  propagating  in  the  circular  waveguide 
were  presented.  In  the  final  report,  we  have  considered  (TEn)  Hu  mode  that  is  the  lowest 
mode  propagating  in  the  circular  waveguide.  The  final  report  is  organized  as  follows.  In 
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Section  2  main  features  of  the  Hu  mode  in  the  circular  waveguide  are  recalled.  The  layout  of 
the  SE  in  the  circular  waveguide,  the  output  signal  and  the  sensitivity  of  the  RS  are 
considered  in  Section  3.  A  finite  difference  time  domain  (FDTD)  method  and  programs  used 
for  electromagnetic  field  components  calculation  in  the  cylindrical  coordinate  system  are 
briefly  presented  in  Section  4.  Modeling  results  for  two  different  configurations  of  the  SE  are 
described  in  Section  5  and  Section  6,  respectively.  In  Appendix,  the  test  results  of  the 
program  for  Hu  mode  are  presented. 

2  Hu  mode  in  a  circular  waveguide 

Features  of  the  Hu  mode  in  the  circular  waveguide  are  summarized  in  this  section. 
Electromagnetic  field  components  as  well  as  the  power  transmitted  through  the  waveguide  are 
considered. 


2.1  Electromagnetic  field  components 

Although  the  H0i  mode,  which  has  been  considered  in  the  interim  report  [3],  is  the 
simplest  TE  mode  in  a  circular  waveguide,  it  is  not  the  lowest  one.  This  illustrates  the  data 
presented  in  Fig.  1,  where  critical  wavelengths  and  frequencies  for  different  modes  in  a 
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mode  is  the  lowest.  The 
regular  Hu  wave  has  five 
components:  E ^  Er ,  H ^  Hr 
and  Hz.  As  one  can  see  from 
Fig.  2,  the  components  are 
dependent  on  the  azimuthal 
angle  cp.  Plane  A-A,  where 
the  electric  field  has  Er 

component  only,  sometimes  is  called  a  polarization  plane. 

Dependences  of  the  amplitudes  of  the  regular  wave  on  the  radial  r  and  azimuthal  cp 
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Critical  frequencies  and  wavelengths  for  different  modes  in 
circular  waveguide. 


Fig.  2.  The  view  of  lines  of  the  electromagnetic  field  in  the  circular  waveguide  for  mode  Hu. 
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(1) 


coordinates  can  be  expressed  in  the  following  way  [4] 

Er(r , (p)  =  Aa  cos cp, 

Mnr 

E  (r,(p)  =  -AJl(/jllr/a)sm(p, 


Hr  (. r,z)  =  -A—  —J,  (/jnr  /  a )  sin  (p , 

4,  VM> 


H  (r,z)  =  Aa  — 

K  v  #> 

/i 


A  js0  Jl (/Jnr /  a) 


Air 


COS^?. 


Hz  (r,  z)  =  A—  4 1— (//ur  /  a)  sin  q> . 
Mo 


(2) 

(3) 

(4) 

(5) 


Here  the  coefficient  A  is  measured  in  electric  field  units  and  it  depends  on  the  power 
propagating  through  the  waveguide,  A,  Aw,  and  Ac  are  the  wavelength  of  electromagnetic 
oscillations  in  free  space,  in  a  waveguide,  and  the  critical  wavelength  for  the  Hu  mode, 
respectively,  So  and  jUq  are  the  vacuum  permittivity  and  permeability,  J\  and  Jx  are  the  first 
kind  and  first  order  Bessel  function  and  its  derivative,  //n=1.8412  is  the  root  of  Jx : 
J\(M\ i)  -  0 ,  and  a  is  a  radius  of  the  waveguide.  The  shape  of  the  Bessel  functions  appearing 
in  (l)-(5)  expressions  is  shown  in  Fig.  3.  It  is  seen  that  the  maximum  of  radial  and  azimuthal 
components  of  the  electromagnetic  field  is  reached  in  the  center  of  the  waveguide  while  the 
axial  component  of  the  magnetic  field  gets  its  maximum  at  a  metal  wall  of  the  circular 
waveguide.  The  following  expression  describes  the  critical  wavelength  for  the  Hu  mode  [4] 


Fig.  3.  The  shape  of  Bessel  functions  appearing  in  (l)-(5). 
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and  the  wavelength  in  the  waveguide  is  expressed  in  a  usual  way 

A.,=  * 


1-(X/XC) 


(7) 


2.2  Power 


The  power  transfer  by  the  Hu  mode  can  by  straightforwardly  computed  by  integrating 
the  Pointing  vector  over  cross  sectional  area  of  the  waveguide.  Making  use  of  (l)-(4)  one  can 
get  the  following  expression  describing  power  transmitted  through  the  waveguide  for  the  Hu 
mode 


p  a2 A2  l£0 

2XwMi  i  V  Mo 


'{  f(x) 

x  J 


M 11  2 

xdx  +  J  (7j(x))  xdx 


(8) 


The  coefficient  A  can  be  related  to  the  maximum  electric  field  amplitude  in  the  regular  Hu 
wave  in  the  center  of  the  waveguide.  From  (1)  and  (2)  follows 

^=2  E0,  (9) 


where  lu>  is  the  maximum  electric  field  strength  of  the  components  /:)>  and  Er  in  the  center  of 
the  waveguide.  Inserting  (7)  and  (9)  into  (8)  and  performing  integration  one  can  get  the 
following  expression 
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L  V 

describing  power  transmitted  through  the  waveguide.  Bessel  function  values  are 
7o(y«i  0=0.3 160  and  7i(//i  0=0.58 19.  From  (10),  the  maximum  electric  field  strength  in  the 
waveguide  at  a  given  power  transmitted  through  the  waveguide  can  be  also  determined. 


3  Resistive  sensor  in  circular  waveguide 

In  this  section,  the  RS  in  a  circular  waveguide  is  considered.  Since  the  possible  layouts 
of  the  SE  are  discussed  in  detail  in  the  interim  report  [3]  in  the  first  subsection  we  present 
only  brief  description  of  the  layout  of  the  plate  with  metal  contacts  used  in  the  calculations. 
The  output  signal  and  sensitivity  of  the  RS  are  considered  in  the  following  two  subsections. 

3.1  Layout  of  the  sensing  element 

As  it  was  already  mentioned  in  the  Introduction,  the  performance  of  the  RS  is  based 
on  the  electron  heating  effect  in  semiconductors.  Thus,  the  SE  is  actually  a  resistor  made  from 
n-type  Si  with  two  ohmic  contacts.  Therefore,  any  plate  inserted  in  the  waveguide  with 
properly  arranged  contacts  could  serve  as  the  SE.  The  cross  sectional  view  of  such  a  plate  in  a 
circular  waveguide  and  its  characteristic  transverse  dimensions  are  shown  in  Fig.  4.  The 
length  of  the  plate  in  the  wave  propagation  direction  is  denoted  as  /.  Depending  on  the 
position  of  the  metal  contact  on  the  plate,  two  configurations  showed  in  Fig.  4  have  been 
proposed.  A  plate  with  the  contacts  on  its  top  and  bottom  was  named  as  a  “vertical” 
configuration  (Fig.  4a),  while  a  plate  with  the  contacts  on  sidewalls  of  the  plate  got  a  name 
“horizontal”  (Fig.  4b).  These  two  names  will  be  used  further  without  quotes  to  denote  the 
particular  configuration  of  the  SE.  The  possible  practical  realization  of  both  configurations 
was  considered  in  the  interim  report  [3]. 
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Fig.  4.  A  plate  of  the  semiconductor  with  metal  contacts  on  the  wall  of  the  circular  waveguide:  a)  contacts  on 
top  and  bottom  surfaces  -  “vertical”  configuration,  b)  contacts  on  the  sidewalls  -  “horizontal” 
configuration. 

It  should  be  noted  that  the  amplitudes  of  electric  field  of  the  regular  wave  depends  on 
azimuthal  angle  cp  (ref.  to  (1)  and  (2)).  At  cp-  0°  the  component  Er  gets  maximum  (1)  while  at 
cp  -  90  -  E(p.  Since  the  average  electric  field  in  the  SE  might  be  dependent  on  the  angle 
between  obstacle  and  polarization  plane,  the  calculations  should  be  performed  for  the  SE 
situated  at  different  angles  cp'm  respect  to  the  polarization  plane,  as  it  is  shown  in  Fig.  4. 

3.2  Output  signal 

An  output  signal  from  the  RS  being  measured  using  a  high  input  resistance  measuring 
unit  can  be  expressed  in  the  following  way 

(11) 

where  Uo  is  a  DC  voltage  drop  on  the  sensor  and  A R/R  is  a  relative  resistance  change  of  the 
SE  in  the  microwave  electric  field.  It  is  well  established  that  in  the  limit  of  the  slightly  heating 
electric  field  the  relative  resistance  change  is  proportional  to  the  square  of  the  electric  field  [2] 

^  =  f(E)2  (12) 

where  0  is  so-called  an  effective  warm-electron  coefficient  and  (E)  is  an  average  amplitude 
of  the  electric  field  in  the  SE.  The  region,  where  expression  (12)  is  valid,  is  called  a  warm- 
electron  region.  In  general,  the  effective  warm  electron  coefficient  is  frequency  dependent  and 
it  decreases  while  the  frequency  increases.  In  the  present  investigation  we  do  not  take  into 
account  this  effect  therefore  our  findings  can  be  used  up  to  and  including  X-band  where  the 
decrease  of  /?  with  frequency  is  less  than  3%. 

Experimental  investigations  have  revealed  that  the  warm-electron  approximation  holds 
well  up  to  electric  field  roughly  1  kV/cm  where  A R/R  is  of  the  order  of  10%.  At  a  higher 
electric  field,  there  is  considerable  deviation  of  the  relative  resistance  change  from  the 
dependence  predicted  by  expression  (12).  It  was  found  that  over  a  wider  range  of  the  electric 
field  strength,  the  resistance  change  is  described  by  the  following  empirical  relation,  with  two 
adjustable  parameters 
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(13) 


Afl_Vl  +  4Qr(E)2-l 

R  2  K 

where  kn  describes  the  deviation  of  A R/R  from  quadratic  dependence  predicted  by  (12). 
Typical  values  of  0  and  kn*  for  n-type  Si  at  a  room  temperature  are  collected  in  Table  1. 


Table  1.  The  parameters  0  and  k,,'  for  n-type  Si  at  T  =  300  K,  when  the  electric  field  is  applied  in  <11 1> 
crystallographic  direction. 


Specific  resistance, 

Q-cm 

0, 

cm2/V2 

k  * 

fy n 

5 

9.0  x  10'8 

3.0 

20 

9.3  x  10"8 

3.4 

200 

10.1  x  10'8 

4.3 

3.3  Sensitivity 


Now  we  consider  a  sensitivity  of  the  RS  in  the  linear  region  where  the  output  signal  is 
proportional  to  the  pulse  power  P  propagating  in  the  waveguide.  Since  the  resistance  change 
of  the  SE  is  a  quantity  indicating  pulse  power  in  the  waveguide,  it  is  convenient  to  define  the 
sensitivity  of  the  RS  as 


A RPR 
P 


(14) 


Definition  (14)  is  not  unique.  Sometimes  the  sensitivity  is  defined  as  a  signal  to  power  ratio. 
Considering  that  the  signal  amplitude  depends  not  only  on  the  DC  voltage  drop  on  the  RS  but 
on  the  input  resistance  of  the  measurement  unit  used  for  the  output  signal  measurement,  the 
proposed  definition  of  the  sensitivity  is  likely  more  acceptable. 

Inserting  (12),  (10),  and  (7)  into  (14)  one  can  get  the  following  expression 

V^o/go  f  (£)V 


(T  = 
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na 


[.-4] 

0(Mn) 

l  Mi  J 

\  -(A/4)2 


(15) 


V  o  y 


describing  the  sensitivity  of  the  RS  in  the  linear  region.  It  is  worthwhile  to  remind  that  Eo  in 
the  obtained  expression  denotes  the  maximum  of  the  electric  field  amplitude  in  the  regular 
Hu  wave  in  the  center  of  the  empty  waveguide.  The  average  electric  field  is  the  only 
unknown  quantity  in  (15).  Thus  determining  it,  the  sensitivity  of  the  RS  in  the  linear  region 
can  be  calculated. 

Except  ( E )  the  only  multiplier  in  (15)  that  influences  the  dependence  of  the  sensitivity 
of  the  RS  on  frequency  is  the  square  root  appearing  in  the  denominator  of  the  expression.  It 
indicates  the  fact  that  even  at  the  same  power  level  transmitted  through  the  waveguide,  the 
electric  field  strength  in  it  decreases  with  frequency  because  of  the  wave  dispersion  in  the 
waveguide.  This  leads  to  the  decrease  of  the  sensitivity. 

Circular  waveguides  with  the  Hu  mode  are  usually  employed  in  the  frequency  range 
A/ Ac  =  0.85-0.6  [5],  From  (15)  it  is  easy  to  get  that  the  sensitivity  in  this  frequency  range 
decreases  by  a  factor  of  1.52.  If  the  average  electric  field  in  the  SE  increases  by  a  factor  1.23 
in  the  same  frequency  range,  it  will  compensate  the  decrease  of  the  sensitivity  due  to  the  wave 
dispersion  in  the  waveguide  and  as  a  result,  the  RS  with  flat  frequency  response  will  be 
developed. 
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Considering  this  fact  and  the  requirements  for  the  RS  formulated  in  the  Introduction, 
we  performed  calculations  of  the  average  electric  field  in  the  SE  with  a  goal  to  find  the  RS 
with  the  optimal  frequency  response.  When  (/:)  is  determined,  the  dependence  of  A R/R  on  P 
can  be  calculated  in  a  wider  dynamical  range  of  the  power  transmitted  through  the  waveguide. 
Inserting  (10)  into  (13),  one  can  get  the  following  expression  describing  the  dependence  of  the 
relative  resistance  change  on  power 


A  R 
R 


2k* 


1  + 


</>'  (ri! 


na 


^o(Ml)  + 


4V 


Mu] 


Jl(Pn) 


\l  Mo  /  £q 


rP~  1 


(16) 


3.4  Conclusions 

Two  configurations  of  the  sensing  elements  of  the  RS  for  the  measurement  of  HPM 
pulses  in  the  circular  waveguide  with  the  Hu  mode  have  been  considered.  Factors  influencing 
the  sensitivity  of  the  RS  in  the  circular  waveguide  were  considered,  and  the  requirements  for 
the  SE  with  flat  frequency  response  were  formulated. 

4  FDTD  method 


The  tremendous  increase  of  the  speed  and  memory  of  personal  computers  allows 
solving  three-dimensional  electromagnetic  problems  at  a  reasonable  time.  We  used  the 
simplest  and  straightforward  finite-difference  time-domain  (FDTD)  method  for  the 
calculation  of  the  electromagnetic  field  components  [6,  7]  in  the  structure  under  investigation. 
In  this  section,  the  model  and  program  used  for  the  calculation  of  the  average  electric  field  in 
the  semiconductor  obstacle  placed  on  the  wall  of  the  circular  waveguide  are  briefly  presented. 


4.1  Maxwell’s  equations 


The  modeled  section  of  the  waveguide  with  the  obstacle  is  shown  in  Fig.  5.  We  used  a 
cylindrical  coordinate  system  and  dimensionless  coordinates  and  time:  r!a,  cp,  zJa,  t-v/a  where 
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Fig.  5.  The  sectional  view  of  the  modeled  circular  waveguide  with  obstacle  in  xOz  plane  (a)  and  x()y  plane  (b). 
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v  is  the  velocity  of  light  in  free  space  and  a  is  a  radius  of  the  waveguide.  In  the  plane  z-Zex^  the 
Hu  type  wave  is  excited.  It  propagates  into  both  sides  from  the  excitation  plane.  The  obstacle 
is  placed  at  one  wavelength  in  the  waveguide  ahead  from  the  excitation  plane  and  at  the  same 
distance  before  the  right  side  of  the  modeled  waveguide  section.  In  the  planes  z— 0  and  z-Zmax 
nonreflecting  boundary  conditions  are  satisfied.  Therefore,  the  waves  traveling  left  from  the 
excitation  plane  as  well  as  reflected  from  the  semiconductor  obstacle  are  absorbed  in  the 
plane  z= 0,  whereas  the  wave  passing  semiconductor  structure  is  absorbed  in  the  plane  z=zmax • 
Due  to  the  reflection  from  the  semiconductor  obstacle  the  partly  standing  wave  is  formed 
between  the  planes  z-Zex  and  z=z0b •  From  the  amplitude  distribution  in  this  area,  the  reflection 
coefficient  was  determined. 

The  regular  Hu  type  wave  has  five  electromagnetic  field  components  (ref.  to  (l)-(5)). 
In  a  vicinity  of  the  SE  the  Ez  component  might  also  appear.  Therefore,  to  determine  the 
average  electric  field  amplitude  in  the  semiconductor  obstacle  Maxwell’s  equations  have  to  be 
solved  computing  all  six  components  of  the  electric  and  magnetic  fields.  Making  use  of 
dimensionless  variables  and  measuring  up  the  magnetic  field  strength  in  electric  field  units 
ZqH,  where  Zo  is  an  impedance  of  free  space,  Maxwell’s  equations  in  the  semiconductor 
obstacle  can  be  written  down  in  a  following  way 

-  =  (VxH-rE)/f,  (17) 

®=-(VxE).  (18) 

where  y—  Z^alp  accounts  for  losses  in  the  semiconductor  structure.  Here  p  and  s  are  the 
specific  resistance  and  relative  dielectric  constant  of  the  semiconductor  obstacle,  it  was 
assumed  that  p  - 1  for  the  entire  simulation  area.  Outside  the  semiconductor  obstacle  y—  0 

and  s-  1 .  Making  use  of  cylindrical  coordinates  for  the  calculation  of  the  curl  of  the  electric 
and  magnetic  field  and  replacing  time  and  coordinate  derivatives  in  (17)  and  (18)  with  finite 
differences,  one  can  obtain  the  set  of  equations  for  updating  new  values  of  the  components 
from  the  older  ones.  This  set  can  be  found  in  the  interim  report  [3]. 

The  grid  of  points  where  the  particular  component  is  computed  is  shifted  at  a  half  of 
step  in  respect  to  each  other  as  it  was  proposed  by  Yee  [6].  Moreover,  electric  and  magnetic 
fields  are  calculated  at  different  time  moments  providing  h 2  accuracy  in  the  calculation  both 
space  and  time  derivatives.  The  details  of  the  application  of  this  technique  to  the  cylindrical 
coordinate  system  can  be  found  in  monograph  [7]. 

The  grid  can  be  chosen  in  such  a  way  that  it  starts  and  finishes  with  the  points  where 
the  electric  field  components  should  be  calculated.  Therefore,  on  the  waveguide  wall  the 
components  E <p  and  Ez  are  zeroed.  Depending  on  the  considered  configuration  of  the  SE  the 
corresponding  components  of  the  electric  field  are  set  to  zero  on  the  metal  contacts:  E ^  and  Ez 
for  the  vertical  configuration,  E,  and  Ez  for  the  horizontal  one.  In  the  planes  z=0  and  z-Zmax 
nonreflecting  boundary  conditions  for  the  components  Er  and  E 9  are  satisfied. 

At  t< 0  there  are  no  electromagnetic  fields  in  the  modeled  section  of  the  waveguide, 
therefore  all  components  of  the  electric  and  magnetic  field  are  set  to  zero.  When  the 
dimensions  of  the  semiconductor  obstacle  is  much  less  than  the  characteristic  dimensions  of 
the  waveguide  its  influence  on  the  wave  propagating  in  the  waveguide  is  comparatively  small. 
In  such  a  case,  filling  the  waveguide  with  the  ordinary  Hu  wave  components  the  stationary 
solution  is  achieved  faster. 

Choosing  the  time  step  the  Courant  criterion  formulated  for  3-D  cylindrical  coordinate 
FDTD  procedure  in  [8]  was  taken  into  account.  Having  in  mind  that  the  components  of  the 
regular  Hu  wave  depends  on  the  azimuthal  angle  <p,  in  general  the  average  electric  field  in  the 
SE  depends  on  the  angle  between  the  polarization  plane  and  SE.  Therefore,  the  dependence  of 
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the  average  electric  field  in  the  SE  on  angle  cp  should  be  considered  as  well  (Fig.  5b).  Let  us 
recall,  that  the  radial  component  of  the  electric  field  of  the  regular  Hu  wave  is  at  its  maximum 
in  the  polarization  plane  xOz  (1)  while  the  azimuthal  -  in  the  perpendicular  plane  yOz  (2). 
When  the  center  of  the  obstacle  is  placed  on  v  or  y- axes  the  distribution  of  electromagnetic 
field  amplitudes  is  symmetrical  in  respect  to  xOz  or  yOz  planes  and  saving  computers  memory 
only  half  of  the  waveguide  window  might  be  modeled. 

When  one  is  dealing  with  FDTD  procedure  in  cylindrical  coordinates  a  numerical 
singularity  is  encountered  at  r  =  0.  To  overcome  the  singularity  we  have  followed  a  method 
proposed  in  [9]  that  is  based  on  the  use  of  the  Cartesian  coordinate  system  in  the  vicinity  of 
r-  0.  Since  the  method  was  described  in  details  in  interim  report  [3]  and  it  is  independent  on 
the  mode  propagating  in  the  waveguide,  we  do  not  consider  it  here. 


4.2  The  program 


The  program  computing  the  electromagnetic  field  components  was  written  using  C++ 
programming  language.  It  works  as  follows.  During  calculations  the  amplitudes  of  the 
particular  component  of  the  electromagnetic  field  at  each  point  of  the  investigated  structure 
was  summed  and  stored  in  additional  arrays.  After  each  period,  the  amplitudes  of  the  electric 
field  components  are  determined  and  the  average  value  in  the  SE  is  calculated.  Obtained  value 
is  compared  with  the  value  calculated  one  period  before.  The  calculations  are  terminated 
when  the  relative  difference  between  these  values  had  been  less  than  the  predetermined  value 
8.  Otherwise,  the  successive  period  has  been  modeled.  The  number  of  periods  that  are 
necessary  to  model  depends  on  the  obstacle  size  and  specific  resistance.  For  small  obstacles, 
the  difference  between  amplitudes  less  than  0.01  can  be  achieved  after  3  periods. 

As  it  was  already  mentioned  in  the  previous  section,  the  electromagnetic  field 
distribution  sometimes  might  be  symmetric  in  azimuthal  direction  in  respect  to  a  plane  laying 
in  the  middle  of  the  obstacle.  In  this  case,  we  used  programs  calculating  the  electromagnetic 
field  amplitudes  in  the  half  of  the  waveguide.  It  is  worth  to  mention  that  computed  results 
have  been  compared  with  those  obtained  using  program  calculating  electromagnetic  field 
components  in  a  whole  waveguide  and  no  difference  between  results  has  been  obtained. 

Calculations  have  been  performed  for  the  waveguide  with  the  inner  radius  a  =  1  cm. 
For  such  a  waveguide  the  critical  wavelength  for  the  Hu  mode  is  Ac  =  3.41  cm  (6)  that 
corresponds  to  the  cutoff  frequency  fc=  8.78  GHz.  We  performed  calculations  starting  from 
fjf-  0.85,  f-  10.3  GHz,  towards  higher  frequencies.  Typical  dimensionless  values  of  steps 
used  in  calculations  were  Az  =  Ar=  0.025,  A  cp  -  1.406°,  At  =  3-10"4,  the  value  of  the  absolute 
error  was  set  8  =  0.01. 

The  written  program  has  been  tested  modeling  the  obstacle  that  totally  fills  the  inner 
gap  of  the  circular  waveguide  window.  Calculated  reflection  coefficient  values  have  been 
compared  with  the  analytical  solution  in  Appendix  1.1.  We  also  compared  the  average  electric 
field  calculated  in  the  obstacle  using  present  program  with  the  results  obtained  using  the 
Cartesian  coordinate  system.  These  results  are  described  in  subsection  6.4.  A  good 
coincidence  between  results  obtained  in  different  ways  proves  the  eligibility  of  the  written 
program  for  the  calculation  of  the  average  electric  field  and  the  reflection  from  the 
semiconductor  obstacle  placed  on  the  wall  of  the  circular  waveguide. 

Since  Maxwell’s  equations  are  linear,  the  solution  obtained  in  the  circular  waveguide 
with  a  radius  a  at  a  frequency /for  the  particular  size  and  specific  resistance  obstacle  (/*,  Z,  p) 
can  be  scaled  to  the  other  frequency  in  the  waveguide  with  radius  a!  by  a  simple  linear 
transformation 


h'  =  h~,  /'  =  /—,  p'  =  p~,  f  =  f-r 
a  a  a  a 


(19) 
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where  primed  symbols  denote  a  new  size,  specific  resistance  and  frequency.  As  it  was  already 
mentioned  in  Subsection  3.2  in  the  present  investigation  we  do  not  take  into  account  the 
dependence  of  0  on  frequency.  This  assumption  is  valid  up  to  and  including  X-band. 
Therefore,  wishing  to  compensate  the  frequency  response  of  the  sensitivity  at  higher 
frequency  bands  the  decrease  of  / ?  with  frequency  should  be  taken  into  account  [10]. 

4.3  Conclusions 

The  model  of  the  semiconductor  obstacle  is  considered  and  the  method  to  solve 
Maxwell’s  equations  in  the  cylindrical  coordinate  system  is  briefly  presented.  The 
performance  of  the  program  written  for  the  calculation  of  the  average  electric  field  in  the 
semiconductor  obstacle  and  the  reflection  coefficient  from  it  is  described. 
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5  Modeling  results:  Vertical  configuration 

In  this  section,  the  modeling  results  for  the  vertical  configuration  (contacts  on  top  and 
bottom  surfaces  Fig.  4a)  are  presented.  In  the  first  subsection,  the  resonances  that  occurring  in 
the  dielectric  SE  inserted  in  the  circular  waveguide  with  Hu  mode  and  the  influence  of 
electrophysical  parameters  of  the  SE  on  the  resonance  position  and  amplitude  are  considered. 
The  dependence  of  the  electric  field  in  the  SE  on  the  position  of  the  polarization  plane  in 
respect  to  the  SE  is  described  in  this  subsection  as  well.  The  resistive  sensor  with  optimized 
frequency  response  is  presented  in  the  second  subsection.  The  dependence  of  the  sensitivity 
on  the  angle  between  the  SE  and  polarization  plane  is  considered  in  the  third  subsection.  The 
peculiarities  of  the  average  electric  field  in  the  SE  as  well  as  in  the  whole  waveguide  section 
containing  a  semiconductor  obstacle  are  considered  in  the  fourth  subsection. 

5.1  Resonances 

The  large  value  of  the  relative  dielectric  constant  of  Si  s-  11.9  enforces  considerable 
decrease  of  the  wavelength  of  the  electromagnetic  wave  in  it.  Therefore,  even  in  a  small 
dimension  obstacles  placed  in  the  waveguide  some  resonances  might  occur.  It  was  also 
confirmed  by  our  earlier  investigations  of  the  semiconductor  obstacles  in  the  rectangular 
waveguide  [10,  11].  In  the  interim  report  [3],  we  have  demonstrated  that  for  Hoi  mode  the 
position  of  the  resonance  in  the  frequency  scale  and  its  type  strongly  depends  on  the 
configuration  of  the  metal  contacts.  Although  the  distribution  of  electromagnetic  field 
components  in  Hu  mode  is  more  complicated,  the  same  peculiarities  as  for  H0i  mode  are 
characteristic  of  the  insulating  obstacle  in  the  waveguide  with  Hu  mode.  Namely,  Er 
component  is  the  largest  in  the  vertical  configuration  SE  while  Eq>  is  the  largest  in  the 
horizontal  one.  The  other  components  of  the  electric  field  in  the  obstacle  are  much  smaller  in 
comparison  with  the  main  in  both  configurations.  Considering  the  distribution  of  the  largest 
electric  field  component  within  obstacle  in  the  plane  qfiz  at  the  resonance  conditions,  it  was 
found  that  the  distribution  is  more  or  less  symmetric  in  respect  to  the  center  of  the  SE  in  the 
wave  propagating  and  transverse  directions.  The  resonance  in  general  is  observed  when  either 
the  length  of  the  structure  Z  or  width  corresponds  to  the  integer  number  of  AJ2,  where  Ad  is  the 
wavelength  of  microwaves  in  the  structure.  Since  the  resonance  might  be  used  as  a  tool  for  the 
engineering  of  the  output  signal  dependence  on  frequency  [10,  11]  we  considered  in  more 
detail  the  influence  of  electrophysical  parameters  of  the  SE  (size  and  specific  resistance)  to 
the  resonance  phenomena  in  it. 

The  survey  of  those  investigations  for  the  vertical  configuration  is  collected  in  Fig.  6 
where  the  dependences  of  the  average  electric  field  in  the  obstacle  on  frequency  for  different 
size  and  specific  resistance  obstacles  are  presented.  The  electric  field  is  normalized  to  the 
maximum  amplitude  of  the  electric  field  in  the  centre  of  the  empty  waveguide.  The  center  of 
the  SE  is  located  in  the  polarization  plane  (cp=  0°,  ref.  to  Fig.  5b)  where  the  component  E,  is 
at  its  maximum  (1).  Initially,  we  investigated  resonances  of  the  dielectric  obstacle  by 
changing  one  of  the  parameters;  say  the  length  Z,  while  the  other  parameters  had  been  fixed. 
The  dependence  of  the  resonance  position  on  Z  is  shown  in  Fig.  6a.  It  is  seen  that  when  the 
length  of  the  SE  is  growing,  the  resonance  shifts  to  the  lower  frequency.  It  is  usual  behavior 
of  the  resonance  since  the  larger  wavelength  wave  fits  the  increased  length  of  the  obstacle. 
The  dependence  of  the  resonance  position  on  the  transverse  dimension  is  shown  in  Fig.  6b.  It 
is  seen  that  the  increase  of  d  also  leads  to  the  decrease  of  the  resonance  frequency.  One  can 
suppose  that  the  increase  of  the  transverse  dimension  of  the  obstacle  leads  to  the  stronger 
coupling  of  it  with  the  wave  propagating  through  the  waveguide  and  this  is  the  reason  shifting 
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Fig.  6.  Calculated  dependences  of  the  average  of  the  component  Er  in  the  vertical  configuration  dielectric 
obstacle  on  frequency  for  the  different  length  /  (a),  width  d  (b),  and  height  h  (c)  obstacles.  The 
dependences  of  (Er)  on  frequency  for  different  specific  resistance  obstacles  are  shown  in  subfigure  (d). 
The  dimensions  of  the  obstacle  and  specific  resistance  (d)  of  it  are  specified  in  the  figure.  Relative 
dielectric  constant  of  the  obstacle  is  set  1 1.9,  a  =  10  mm. 

the  resonance  to  the  lower  frequency.  The  results  of  the  calculations  for  the  different  height 
obstacles  are  shown  in  Fig.  6c.  It  is  seen  that  the  increase  of  h  leads  to  the  shift  of  the 
resonance  toward  higher  frequency.  It  is  worth  mentioning  that  this  shift  is  not  so  large. 
Changing  the  height  of  the  obstacle  three  times,  the  resonance  is  shifted  only  by  0.9  GHz.  In 
the  plane  transverse  to  the  axial  axis,  the  obstacle  looks  as  the  sector  of  the  ring.  While  the 
height  of  the  SE  increases,  the  actual  width  of  it  decreases.  One  might  expect  that  this 
decrease  of  the  actual  width  is  responsible  for  the  shift  of  the  resonance  to  the  higher 
frequency.  Introduction  of  losses  in  the  system  usually  dumps  resonance.  The  SE  in  the 
circular  waveguide  is  not  the  exception.  As  follows  from  the  results  presented  in  Fig.  6d, 
electric  field  strength  considerably  decreases  when  the  conductivity  of  the  SE  increases. 

The  angular  dependence  of  the  resonance  in  the  SE  was  also  investigated  since  the 
amplitudes  of  the  electric  field  in  the  regular  Hu  wave  depends  on  cp.  When  cp-  0°  (ref.  to 
Fig.  5b)  the  obstacle  is  mainly  exposed  by  Er  component  since  E9  =  0  (2).  At  cp  -  90°,  on  the 
contrary,  the  SE  is  mainly  influenced  by  E 9  since  Z?r  =  0  (1). 

Calculated  dependences  of  the  averaged  electric  field  in  the  SE  on  frequency  for 
different  cp  are  shown  in  Fig.  7.  As  one  can  see  from  the  figure,  for  the  strip  shape  obstacle 
(Fig.  7a)  the  position  of  the  resonance  is  independent  of  the  azimuthal  angle.  The  average 
electric  field  component  in  the  SE  decreases  while  the  obstacle  is  shifted  towards  the  y-axis. 
Considering  distribution  of  the  component  Er  inside  SE  in  cpttz  plane  it  was  found  that  the 
resonance  appears  with  the  maximal  electric  field  at  the  beginning  and  end  of  the  obstacle  in  z 
direction  and  minimal  electric  field  in  the  midpoint  of  it.  It  seems  that  such  type  of  the 
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Fig.  7.  Calculated  dependences  of  the  average  of  the  component  Er  in  the  dielectric  obstacle  on  frequency  for 
the  strip  shape  (a)  and  plate  shape  type  (b)  obstacles  located  at  angle  cp  from  the  x-axis.  The  dimensions 
of  the  obstacle  and  shift  angle  are  specified  in  the  figure.  Relative  dielectric  constant  of  the  obstacle  is 
set  11.9,  a  =  10  mm. 
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Fig.  8. 


resonance  is  tightly  connected  with  the  Er 
component  of  the  regular  wave.  Moving  the 
SE  towards  y-axis  the  amplitude  of  the  Er  in 
regular  wave  decreases  causing  in  turn  the 
decrease  of  the  average  electric  field 
amplitude  in  the  obstacle. 

The  dependences  of  the  (Er)  on 
frequency  are  more  complicated  for  the  plate 
shape  obstacle  which  dimensions  in 
longitudinal  ant  transverse  direction  are 
comparable.  As  it  is  seen  from  Fig.  7b,  at  a 
cp-  0  the  resonance  appears  at  11.0  GHz. 

While  the  obstacle  is  shifted  towards  the  y- 
axis,  this  resonance  disappears  and  a  new 
one  emerges  at  a  frequency  12.4  GHz.  From 
the  distribution  of  the  electric  field 
component  Er  in  a  cp  0 z  plane  of  the  obstacle 
shown  in  Fig.  8,  it  is  seen  that  the 

competition  between  longitudinal  and  transverse  resonances  takes  place  in  the  plate  shape 
obstacle.  When  the  obstacle  is  mainly  exposed  by  the  component  Er  (<p=  0),  the  longitudinal 
resonance  occurs.  It  seems  that  transverse  resonance  appears  when  the  obstacle  is  exposed  by 
the  component  E ^  At  some  intermediate  angle,  both  resonances  can  be  distinguished  in  the 
dependence  of  the  (Er)  on  frequency  (Fig.  7b). 

In  conclusion,  it  may  be  said  that  performed  investigations  revealed  the  ways  to 
manage  the  resonance  that  occurs  in  the  SE  placed  on  the  wall  of  the  circular  waveguide  with 
Hu  mode.  As  in  a  case  of  H0i  mode  the  position  of  the  resonance  in  the  frequency  scale  can 
be  easily  shifted  by  changing  either  the  length  or  the  width  of  the  obstacle  while  the  average 
electric  field  in  it  can  be  adjusted  by  changing  specific  resistance  of  the  SE.  Hu  mode  is  more 
complicated  since  the  resonance  might  depend  on  the  position  of  the  SE  in  respect  to  the 
polarization  plane.  Therefore,  in  our  further  investigations  for  searching  the  optimal  frequency 
response  we  shall  start  from  the  SE  located  symmetrically  in  respect  to  the  polarization  plane. 
Angular  dependence  of  the  sensitivity  will  be  considered  afterwards. 


Distributions  in  a  cp  0z  plane  of  the  component 
Er  in  the  obstacle  at  a  resonance  conditions: 
cp=0°,  f=  11.0  GHz,  cp=  45°,  /=  12.4  GHz. 
Dimensions  of  the  obstacle  are  indicted  in  the 
figure. 
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a)  h  -  1  mm,  p-  10  ncm 
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Fig.  9.  The  ratio  of  the  average  electric  field  of  the  component  Er  calculated  at  14.5  GHz  to  the  (Er)  calculated 
at /=  10.3  GHz  versus  longitudinal  and  transverse  dimensions  of  the  SE  for  p  =  10  Q-cm  (a)  and 
p=  20  Q-cm  (b).  The  number  near  the  point  corresponds  to  the  value  of  the  ratio  for  the  particular  pair 
of  d  and  /.  Waveguide  radius  is  10  mm,  the  height  of  the  SE  is  set  1  mm.  The  solid  curve  in  the  figures 
corresponds  to  d  and  /  at  which  the  DC  resistance  of  the  RS  is  50  Q.  Therefore,  above  this  curve,  the 
resistance  of  the  RS  is  less  than  50  Q;  below  this  curve,  it  is  more  than  50  Q. 


5.2  Frequency  response  optimization 

To  optimize  the  frequency  response  of  the  RS,  we  calculated  the  dependence  of  the 
average  electric  field  < Er )  on  frequency  in  the  range  10.3-14.5  GHz  that  corresponds  to  A// lc  = 
0.85-0.6.  This  frequency  range  has  been  considered  in  Subsection  3.3  and  it  was  shown  that 
the  increase  of  (Er)  by  a  factor  1.23  is  necessary  compensating  the  decrease  of  the  sensitivity 
due  to  wave  dispersion  in  the  waveguide.  Therefore,  our  task  is  to  find  the  values  of  h,  d,  l, 
and  p  at  which  such  increase  of  (Er)  with  frequency  is  observed.  As  it  was  already  mentioned 
in  the  previous  subsection  initially  we  investigated  the  SE  placed  on  the  x-axis  where  Er  is  at 
its  maximum. 

Taking  into  account  the  requirements  for  the  RS  formulated  in  Introduction,  we  started 

our  search  for  the  set  of  the  optimal  electrophysical  parameters  of  the  SE  considering  its 

height  h  =  1  mm.  Such  height  comprises  1/10  of  the  radius  of  the  waveguide  so  one  might 

expect  small  perturbation  of  the  waveguide  by  the  obstacle.  Results  of  the  calculations  are 

shown  in  Fig.  9.  Numbers  in  the  figure 

near  the  particular  point  show  a  ratio 

(£V)(14.5  GHz)/<£r>(10.3  GHz)  calculated 

for  the  particular  pair  of  d  and  /.  Solid  line 

in  the  figures  corresponds  to  50  Q  DC 

resistance  of  the  RS.  Below  this  curve,  the 

resistance  of  the  RS  is  larger  than  50  Q, 

whereas  above  it  R  <  50  Q.  Considering 

results  shown  in  the  figure,  one  can  see 

that  in  general  much  larger  increase  of  ( Er ) 

with  frequency  is  observed  than  it  is 

necessary  to  compensate  decrease  of  the 

electric  field  in  the  waveguide  due  to  wave 

c.  in  „  .  f/rvucru  dispersion.  Only  very  short  sensors 

Fig.  10.  Dependence  of  (£r)(14.5  GHz)/(£r)(10.3  GHz)  on  1  \  J 

longitudinal  and  transverse  dimensions  of  the  SE  demonstrate  the  increase  of  ( Er )  that  is 

for  p=2Q  Q-cm,  h  =  0.5  mm,  a  =  \  cm.  more  or  less  acceptable  for  the 


b  -  0.5  mm,  p  -  20  Qcm 
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compensation  of  the  frequency  response,  but  they  resistance  is  larger  than  50  Q  and  such  a 
sensor  is  not  suitable  for  the  registration  of  the  short  HPM  pulses.  It  should  be  noted  that  the 
situation  is  not  improved  by  decreasing  specific  resistance  of  SE.  Calculations  for  p  =  5  Q-cm 
demonstrate  practically  the  same  growth  of  ( Er )  within  considered  frequency  range.  It  seems 
that  the  SE  with  p-  20  Q-cm  and  Z  =  5  mm  demonstrate  the  smallest  increase  of  electric  field 
with  frequency  but  it  resistance  is  more  than  50  Q. 

The  resistance  of  the  sensor  decreases  by  decreasing  its  the  height.  Thus,  the  SE  with 
h  -  0.5  mm  was  investigated  expecting  to  get  the  sensor  with  desirable  increase  of  the  average 
electric  field  with  frequency  and  resistance  lower  than  50  Q.  Calculation  results  for  such  SE  is 
shown  in  Fig.  10.  As  one  can  see  from 
the  figure,  for  each  width  of  the  obstacle 
except  the  widest  one  (d  =  4.9  mm)  in 
the  length  range  4. 5-5. 5  mm  at  least  a 
few  points  can  be  found  where  the 
calculated  ratio  of  (Er)(14.5 
GHz)/(£V)(10.3  GHz)  is  close  to  the 
desirable  value  1.23.  Detailed 
investigations  of  the  frequency  response 
for  different  width  d  and  length  Z  of  the 
SE  were  performed.  The  SE  under 
investigation  was  located  in  the 
polarization  plane,  hence  <p=  0.  The  results  of  calculation  are  collected  in  Table  2.  In  the  first 
column  of  the  table  the  width  of  the  SE  is  presented.  Calculations  for  particular  d  have  been 
performed  searching  for  Z  that  provides  minimum  sensitivity  variation.  The  length  and 
minimum  sensitivity  variation  are  displayed  in  the  second  and  third  columns  of  the  table. 

As  one  can  see  from  the  table,  the  smallest  variation  of  the  sensitivity  is  characteristic 
of  the  SE  with  dimensions  d  =  3.8  mm  and  Z  =  5.0  mm.  Its  DC  resistance  is  roughly  21  Q.  At 
a  moment  it  is  the  optimal  SE  demonstrating  the  smallest  variation  of  the  sensitivity  at  <p  =  0. 
The  dependence  of  the  sensitivity  on  frequency  calculated  with  the  help  of  (15)  for  such  RS  is 
shown  in  Fig.  1  la.  In  the  calculation  the  value  of  0  listed  in  Table  1  was  used.  As  one  can  see 
from  the  figure,  flat  frequency  response  is  characteristic  of  the  optimal  RS.  Due  to  a  small 
height,  the  optimal  sensor  does  not  perturb  much  the  field  distribution  in  the  waveguide. 


Table  2.  The  length  /  of  the  SE  at  which  calculated  sensitivity 
variation  is  at  its  minimum  within  10.3-14.5  GHz 
frequency  range  for  different  width  d :  h  =0.5  mm, 
p  —  20  Q-cm,  a  =  1  cm,  cp  =  0°. 


d,  mm 

Z,  mm 

A  C'C,  % 

1.9 

5.5 

±7.0 

2.4 

5.5 

±6.0 

2.9 

5.0 

±6.0 

3.4 

5.0 

±4.5 

3.8 

5.0 

±3.5 
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Fig.  11.  Calculated  dependence  of  the  sensitivity  on  frequency  (a)  and  dependence  of  the  relative  resistance 
change  on  power  transmitted  through  the  waveguide  at/=  12  GH  for  the  optimal  RS.  Solid  line  in  (b) 
shows  results  calculated  using  (16)  and  dotted  line  demonstrates  linear  dependence  that  is  characteristic 
of  the  warm  electron  region.  Dimensions  of  the  SE  are  hxdxl  =  0.5x3. 8x5  mm3, 
p  —  20  Q-cm,  a  -  10  mm,  cp  =  0°,  R  =  21  Q. 
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hxdxl  =  0  5x3  8x5  mmJ,  p  =  20  Dcm 


Fig.  12.  Dependences  of  the  sensitivity  on  frequency  for 
the  optimal  SE  at  arbitrary  position  in  respect  to 
the  polarization  plane.  An  angle  between  it  and  the 


Calculated  VSWR  slightly  increases  with 
frequency,  but  it  is  less  than  1.03  within 
considered  frequency  range. 

Once  the  electric  field  strength  in 
the  SE  is  determined,  the  dependence  of 
the  relative  resistance  change  on  a  power 
transmitted  through  the  waveguide  was 
calculated  using  expression  (16)  and 
parameters  of  n-Si  presented  in  Table  1. 
First,  we  determined  maximum  power  up 
to  which  the  waveguide  with  a  radius 
a  -  10  mm  can  be  potentially  used. 
Considering  the  electric  field  strength 
25000  V/cm  as  a  limiting  due  to  a 
breakdown  of  air,  one  can  calculate  the 
transmitted  maximum  power  Pmax  for  the 
Hu  mode  using  expression  (10).  At 


center  of  the  SE  is  indicated  in  the  figure:  f—  10.3  GHz  0.65  MW  was  obtained.  As 
hxdxl  =  0.5x3. 8x5  mm3,  p  =  20  Q-cm,  <2=10  mm,  it  was  already  mentioned  (Subsection 
^  =  3.3),  at  higher  frequency  more  power  is 

required  to  get  the  same  electric  field 
strength.  Therefore,  one  can  consider  that  the  maximum  power  that  can  be  transmitted 
through  the  waveguide  (a  =  10  mm)  without  breakdown  is  roughly  0.65  MW  and  the 
dependence  of  the  relative  resistance  change  on  pulse  power  transmitted  trough  the 
waveguide  is  calculated  using  expression  (16).  Calculation  results  are  shown  in  Fig.  lib  by  a 
solid  line.  Dotted  line  corresponds  to  the  linear  dependence  of  A R/R  on  P  that  is  characteristic 
of  the  warm  electron  region.  It  is  seen  that  at  a  maximum  power  roughly  40%  relative 
resistance  change  of  the  SE  is  observed.  It  is  worthwhile  mentioning  that  the  change  of 
resistance  in  general  should  influence  the  average  electric  field  in  the  SE.  Therefore,  the  curve 
in  Fig.  lib  in  a  high  power  limit  should  be  considered  as  some  approximation  useful  to 
estimate  relative  resistance  change  of  the  RS  in  the  HPM  region. 

It  should  be  noted,  that  the  optimal  RS  could  be  used  in  the  polarization  plane  only. 
When  the  SE  is  shifted  out  of  plane,  its  sensitivity  becomes  dependent  on  frequency.  It 

illustrates  calculation  results  shown  in 
Fig.  12.  It  is  seen  that  by  shifting  the  SE 
from  the  polarization  plane  its  sensitivity 
increases.  The  increase  of  the  sensitivity 
strongly  depends  on  frequency. 
Comparing  sensitivity  at  <p=0°  and 
cp  -  90°  it  is  seen  that  at  a  lower 
frequency  the  sensitivity  increases  by  a 
factor  6.5  while  in  a  high  frequency  range 
-  by  1.9  only.  As  a  result,  at  a  fixed  angle 
cp  the  considerable  decrease  of  the 
sensitivity  with  frequency  is  observed. 

Either  polarization  plane  is 
unknown,  or  it  can  be  turned  off  by  the 
accidental  deformation  of  the  waveguide 
[4],  the  measurement  accuracy  employing 
the  optimal  RS  might  be  poor.  One  of  the 


h  -  0.5  mm,  p~  20  Ocm,  90* 


Fig.  13. 


f,  GHz 

Dependences  of  the  sensitivity  on  frequency  for 
the  different  cross-sectional  SE  placed  at  cp=  90°. 
Cross-sectional  size  of  the  SE  is  indicated  in  the 
figure:  h  -  0.5  mm,  p=  20  Q-cm,  a=  10  mm. 
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3 )  hxdxl  -  0.5x2  9x5  mm3,  p-  20  Qcm 


Fig.  14.  Calculated  dependences  of  the  sensitivity  on  frequency  (a)  for  several  cp  and  dependences  of  the 
sensitivity  on  cp  for  several  frequencies.  The  angles  and  frequencies  are  indicated  in  the  figure.  Solid 
lines  in  (b)  corresponds  to  the  third  order  polynomial  approximation.  Dimensions  of  the  SE  are 
hxdxl  =  0.5x2.9x5  mm3,  p-  20  Q-cm,  a  =  10  mm,  R  =  28  Q. 

possible  solutions  of  the  problem  is  the  design  of  the  SE  that  can  be  rotated  in  the  plane 
perpendicular  to  the  axial  direction.  By  turning  the  SE,  the  position  can  be  found  where  the 
output  signal  is  at  its  minimum.  This  position  corresponds  to  the  polarization  plane  in  the 
cylindrical  waveguide.  Therefore  using  such  design  polarization  plane  can  be  found 
experimentally.  In  the  next  subsection,  the  other  possibility  to  solve  the  problem  using  more 
than  one  SE  in  the  waveguide  will  be  considered. 


5.3  Angular  dependence  compensation 

Considering  angular  dependences  of  the  sensitivity  of  the  SE  listed  in  Table  2,  it  was 
found  that  decreasing  its  width  flatness  of  the  sensitivity  is  improved.  It  illustrates  calculation 
results  shown  in  Fig.  13  for  the  SE  situated  at  <p=  90°  from  the  polarization  plane.  As  one  can 
see  from  the  figure,  a  steep  decrease  of  the  sensitivity,  that  is  characteristic  of  the  widest  SE 
under  investigation,  changes  to  the  increase  while  the  width  of  the  SE  is  reduced.  For  the 
further  investigation  we  choose  the 
SE  with  dimensions  dxl  =  2.9x5 
mm  demonstrating  the  smallest 
±10%  sensitivity  variation  within 
considered  frequency  range. 

We  calculated  the 
dependences  of  the  sensitivity  on 
frequency  for  such  size  SE  for 
different  angle  cp  between 
polarization  plane  and  center  of  the 
SE.  These  calculations  results  are 
shown  in  Fig.  14.  As  one  can  see 
from  the  figure,  when  the  angle  (p 
increases,  the  sensitivity  charac¬ 
teristic  becomes  overcompensated  - 
electric  field  strength  in  the  SE 
increases  more  than  it  is  necessary 
to  compensate  the  wave  dispersion 
in  the  waveguide.  The  dependences 
of  the  sensitivity  on  the  angle 


Fig.  15.  Dependences  of  the  summed  sensitivity  of  two  SE  on  the 
angle  between  one  of  the  SE  and  polarization  plane  for 
several  frequencies.  Lines  correspond  to  the  results 
obtained  using  polynomial  approximation,  points  -  direct 
calculation.  In  the  inset,  the  layout  of  the  sensors  in  the 
waveguide  is  shown  schematically.  Parameters  of  the  SE: 
hxdxl  -  0.5x2.9x5  mm3,  p  -  20  Q-cm,  R  =  28  Q. 
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hxdxl=d5x3,8xS  mm,  f=1 4,6  GHz 


a) 


hxdxb0.5x3.8x5  mm3,  f=10,3  GHz 


b) 


E  IE, 


Fig.  16.  Calculated  distribution  of  the  averaged  in  r  direction  electric  field  component  Er  within  SE  at  10.3  GHz 
(a)  and  14.5  GHz  (b).  Dimensions  and  specific  resistance  of  the  SE  are:  hxdxl  =  0.5x3. 8x5  mm3, 
p  -  20  Q-cm,  it  is  placed  at  cp  =  0°,  a  -  10  mm. 

between  the  center  of  the  SE  and  polarization  plane  are  shown  in  Fig.  14b.  Points  in  the  figure 
correspond  to  the  calculation  results;  lines  demonstrate  approximation  of  these  results  by  third 
order  polynomial.  It  is  seen  that  the  approximation  fits  well  the  calculation  results. 
Considering  in  more  detail  the  dependences  of  the  sensitivity  on  cp  one  can  see  that  they  are 
really  symmetric  in  respect  to  cp=  45°.  This  fact  suggests  an  idea  to  employ  two  SE  placed  at 
cp  -  90°  in  between  and,  by  summing  up  their  signals,  compensate  the  angular  dependence  of 
the  sensitivity.  The  idea  was  examined  using  analytical  expressions  of  £(cp)  obtained 
approximating  calculation  results  by  third  order  polynomial  and  having  in  mind  that  by 
changing  cp'm  the  range  0-90°  the  sensitivity  of  the  SE  number  2  shown  in  the  inset  in  Fig.  15 
will  be  4(90 °-cp).  Calculated  in  such  a  way  dependences  of  the  summed  sensitivity  on  cp  for 
several  frequencies  are  shown  in  Fig.  15  by  lines.  Results  of  the  direct  calculation  of  two  SE 
in  the  waveguide  are  shown  in  the  same  figure  by  points.  It  is  seen  that  making  use  of  two 
sensors  the  angular  characteristic  of  the  RS  is  significantly  improved.  Calculation  shows  that 
using  two  sensors  placed  at  90°  from  each  other  and  measuring  the  sum  of  their  signals,  one 
can  obtain  readings  independent  of  the  position  of  the  polarization  plane  in  respect  to  the  SE 
with  the  accuracy  of  ±7%.  It  is  worthwhile  mentioning  that  the  reflection  from  two  sensors 
increases  negligibly  in  comparison  with  single  one.  Calculated  value  of  the  VSWR  was  less 
than  1.05  in  the  considered  frequency  range. 

5.4  Electric  field  in  the  sensing  element 

To  have  an  idea  about  the  electric  field  strength  in  the  SE  the  field  distribution  of  the 
component  Er  is  shown  in  Fig.  16  for  the  SE  with  the  smallest  sensitivity  variation  found.  The 
dimensions  of  the  obstacle  are  hxdxl  =  0.5x3. 8x5  mm  ,  its  specific  resistance  p=  20  Q-cm, 
and  it  is  placed  at  cp  -  0°.  Er  component  is  averaged  within  obstacle  in  r  direction.  It  is  seen 
that  Er  distribution  is  symmetric  in  respect  to  the  azimuthal  direction.  For  the  lower 
frequency,  the  maximum  of  electric  field  is  formed  at  the  beginning  of  the  SE  and  it  shifts  to 
the  centre  when  frequency  increases.  The  same  dependence  is  characteristic  of  the  narrower 
SE  used  in  the  two  SE  configuration.  It  confirms  results  shown  in  Fig.  17a  and  b  where  the 
distribution  of  Er  at  cp  -  0°  is  shown.  At  cp  -  45°,  the  distribution  of  Er  in  the  SE  becomes 
asymmetric  forming  electric  field  maximum  in  the  opposite  corners  of  the  structure  (ref.  to 
Fig.  17c  and  d).  When  the  SE  is  situated  on  y-axis  (cp=  90°),  it  is  seen  from  Fig.  17e  and  f  that 
distribution  of  the  component  Er  is  again  symmetric  with  a  minimum  in  the  center  of  the  SE 
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c)  E  /E  hxdxl=0.5x2.9x5  mm3,  f=10.3  GHZ.  45& 

r  0 


d) 


E  /E 
r  0 


hxdxl=Q.Sx2.9xS  mm3,  f=14.5  GHz,  45° 


e) 


hxdx  1=0.5x29x5  mm3,  f=10.3  GHz,  90° 


Fig.  17.  Calculated  distributions  of  the  averaged  in  r  direction  electric  field  component  Er  within  SE  at 

10.3  GHz  (a),(c),(e)  and  at  14.5  GHz  (b),(d),(f);  (p=  0°  (a),(b),  cp=A5°  (c),(d),  cp=90°  (e),(f). 

Dimensions  and  specific  resistance  of  the  SE  are:  hxdxl  =  0.5x2. 9x5  mm3,  p=  20  Q-cm,  a  =  10  mm. 

in  the  azymuthal  direction.  The  change  of  the  symmetry  in  Er  distribution  in  the  obstacle  was 
mentioned  in  subsection  5.1  considering  distribution  of  the  electric  field  in  the  dielectric 
obstacle  under  resonance  conditions. 

Electric  field  in  the  SE  mainly  has  the  Er  component.  The  other  two  components  E(p 
and  Ez  are  much  smaller.  It  illustrates  calculation  results  of  components  of  the  average  electric 
field  in  the  SE  for  10.3  and  14.5  GHz  and  several  angles  of  cp  shown  in  Table  3  for  the  SE 
with  dimensions  and  specific  resistance:  hxdxl  =  0.5x2. 9x5  mm3,  p-  20  Q-cm.  As  follows 
from  the  table,  ( E ^  is  roughly  200  times  less  than  (Er)  while  (Ez)  is  even  smaller.  It 
worthwhile  mentioning  that  the  increase  of  h  leads  to  some  grow  of  E ^  and  Ez. 


21 


Table  3.  The  average  electric  field  components  in  the  SE  for  two  frequencies  and  several  azimuthal  angles  cp. 
Dimensions  of  the  SE  hxdxl  =  0.5x2.9x5  mm2,  p  =  20  Q-cm. 


<p,° 

10.3  GHz 

14.5  GHz 

(Er)/Eo 

(E9)/Eo 

(Ez)/Eo 

( Er)/E0 

(E0)/Eo 

<- Ez)/E0 

0 

0.1401 

0.0006 

0.0004 

0.1641 

0.0008 

0.0006 

45 

0.1950 

0.0010 

0.0003 

0.2568 

0.0013 

0.0005 

90 

0.2339 

0.0013 

0.0002 

0.3131 

0.0016 

0.0003 

z/a 
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Fig.  18.  Distributions  of  components  of  the  electric  field  in  the  waveguide  with  the  SE  on  the  different  surfaces 
(jpQz  at  f=  10.3  GHz:  at  cp  =  0°  (a)  and  cp  =  90°  (d)  for  r  =  (/-l)Ar  ;  Er  at  cp  =  0°  (b)  for  r  =  (z-0.5)Ar; 
Ez  at  <p  —  0°  (c)  for  r=(z-0.5)Ar;  Calculation  parameters:  Ar  =0.025,  Az=0.05,  A<p=2.8°, 
At  =  4.50- 10"4.  Dimensions  of  the  SE  are  hxdxl  =  0.5x2.9x5  mm3,  p=  20  Q-cm,  a  =  10  mm. 

Distributions  of  components  of  the  electric  field  in  the  waveguide  section  with  the  SE 
are  shown  in  Fig.  18.  Distributions  of  the  components  Eq,  (a)  and  (d),  Er  (b)  and  Ez  (c)  in 
different  (pOz  planes  are  presented.  It  should  be  noted  that  at  cp  -  0°  the  SE  is  actually  situated 
at  cp-  180°  and  cp-  90°  corresponds  to  the  angle  270°.  Those  positions  are  symmetric  in 
respect  to  the  distribution  of  the  electric  field  in  the  regular  wave.  As  it  is  seen  from  the 
figure,  due  to  a  small  height  of  the  SE,  the  perturbation  of  the  regular  components  of  Hu 
wave  is  small  except  the  planes  close  to  the  SE.  The  component  Ez  that  is  absent  in  the  regular 
wave  is  excited  but  its  amplitude  is  very  small;  less  than  1%  in  the  center  of  the  waveguide.  It 
becomes  more  significant  in  the  obstacle  or  near  it.  It  should  be  noted  that  the  distribution  of 
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the  electric  field  does  not  changes  much  with  frequency.  Therefore,  the  distributions  only  at 
the  lowest  considered  frequency  are  shown. 

5.5  Conclusions 

The  resonance  occurring  in  the  dielectric  obstacle  with  metal  contacts  on  its  top  and 
bottom  (vertical  configuration)  placed  in  the  circular  waveguide  with  regular  Hu  wave  was 
considered.  It  was  found  that  the  resonance  position  in  a  frequency  scale  could  be  shifted  to 
the  lower  frequency  by  increasing  the  length  and  width  of  the  obstacle.  The  average  electric 
field  strength  in  the  obstacle  can  be  reduced  by  increasing  its  conductivity. 

It  was  established  that  the  average  electric  field  in  the  SE  mainly  consists  of  the  Er 
component.  The  components  Ep  and  E-  are  less  significant.  They  increase  when  the  height  of 
the  SE  grows. 

The  optimal  SE  found  at  a  moment  should  have  the  following  dimensions  and  specific 
resistance:  hxdxl  -  0.5x3. 8x5  mm3,  p-  20  Q-cm.  It  should  be  placed  in  the  polarization  plane 
of  Hu  wave.  Calculated  sensitivity  variation  of  the  optimal  RS  was  within  ±3.5%  in  the 
frequency  range  10.3-14.5  GHz.  The  DC  resistance  of  the  optimal  RS  was  21  Q.  Reflection 
introduced  in  the  waveguide  by  the  RS  is  sufficiently  small:  VSWR  <  1.03  within  frequency 
range  under  consideration.  The  optimal  RS  can  be  used  for  the  measurement  of  HPM  pulses 
propagating  in  the  circular  waveguide  close  to  the  breakdown  in  the  waveguide. 

Since  the  signal  from  the  optimal  RS  depends  on  the  position  of  it  in  respect  to  the 
polarization  plane,  the  SE  rotatable  in  the  plane  perpendicular  to  the  axial  direction  has  been 
proposed.  The  minimum  signal  from  the  RS  corresponds  to  the  position  of  the  polarization 
plane  in  the  waveguide. 

The  alternative  RS  has  been  proposed  consisting  of  two  SE  ( hxdxl  =  0.5x2.9x5  mm3, 
p-  20  Q-cm)  placed  at  90°  from  each  other.  It  was  shown  that  the  sum  of  their  signals  within 
±7%  is  independent  of  the  position  of  the  polarization  plane  in  the  waveguide.  The  DC 
resistance  of  each  sensor  was  28  Q.  Calculated  value  of  the  VSWR  for  two  SE  was  less  than 
1 .05  within  the  considered  frequency  range. 
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6  Modeling  results:  Horizontal  configuration 

In  this  section,  the  modeling  results  for  the  horizontal  configuration  (contacts  on  the 
sidewalls  as  shown  in  Fig.  4b)  are  presented.  In  the  first  subsection,  the  influence  of 
electrophysical  parameters  of  the  SE  on  the  resonance  position  is  considered.  The  influence  of 
specific  resistance  and  dimensions  on  the  frequency  response  of  the  RS  is  considered  in  the 
second  subsection.  The  resistive  sensor  with  optimized  frequency  response  is  described  in  the 
third  subsection.  In  the  fourth  subsection,  the  results  of  modeling  of  the  SE  in  Cartesian 
coordinate  system  are  presented. 

6.1  Resonances 

In  the  interim  report  [3]  considering  the  horizontal  configuration  obstacle  (contacts  on 
the  sidewalls)  in  the  cylindrical  waveguide  with  Hoi  mode,  it  was  shown  that  the  major 
electric  field  component  in  the  SE  is  E ^  The  same  peculiarity  is  characteristic  of  the  obstacle 
in  the  cylindrical  waveguide  with  Hu  mode.  Since  the  maximum  electric  field  in  the  SE  is 
found  at  <p=  90°  (ref.  to  Fig.  5b)  where  E '9  component  of  the  regular  wave  is  at  its  maximum 
(2),  all  further  calculations  presented  here  for  the  horizontal  configuration  are  at  cp  -  90°, 
unless  different  angle  is  indicated. 

For  the  obstacle  placed  at  cp-  90°  dependencies  of  the  resonance  position  on  the 
electrophysical  parameters  of  it  are  the  same  as  for  the  obstacle  in  the  waveguide  with  H0i 
mode  [3].  As  usual,  the  increase  of  the  length  leads  to  the  shift  of  the  resonance  position  to  the 
lower  frequency  and  the  increase  of  the  specific  conductivity  -  to  the  decrease  of  the  electric 
field  strength  in  the  SE.  The  position  of  the  resonance  might  also  be  changed  by  varying  the 
height  of  the  obstacle.  Increase  of  h  leads  to  the  shift  of  the  resonance  to  the  lower  frequency 
while  the  increase  of  the  width  does  not  influence  much  the  position  of  the  resonance. 

We  performed  additional  simulations  to  elucidate  the  influence  of  the  position  of  the 
polarization  plane  of  the  incident  wave  in  respect  to  the  center  of  the  SE  on  resonance 
position  and  amplitude.  Results  of  the  simulation  are  shown  in  Fig.  19.  It  is  seen  that  a  clear 
resonance  was  obtained  at  13.1  GHz.  The  amplitude  of  the  electric  field  decreases  when  the 
polarization  plane  approaches  to  the  SE.  Averaged  in  the  azimuthal  direction  distribution  of 
the  electric  field  at  cp-  90°  is  shown  in  Fig.  19b.  It  is  seen  that  the  AJ2  resonance  in  the  wave 
propagation  direction  is  excited  in  the  structure. 


^  hxdxl  =  4x0.39x4  mm3,  p  =  oo 


f,  GHz 


Fig.  19.  The  average  electric  field  in  the  dielectric  obstacle  versus  frequency  for  several  polarization  angles  of 
the  incident  wave  (a)  and  the  distribution  of  the  electric  field  averaged  in  cp  direction  within  the  SE  at 
the  resonance  frequency  13.1  GHz  and  cp  =  90°.  Dimensions  of  the  obstacle  are  indicated  in  the  figure. 
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It  should  be  noted  that  the  similar  size  obstacle  was  also  investigated  in  vertical 
configuration  (Fig.  7b)  and  two  resonances  were  discovered.  However,  in  horizontal 
configuration  material  of  the  device  is  surrounded  by  three  metal  interfaces  (two  contacts  plus 
a  wall  of  the  waveguide)  instead  of  two  in  vertical  one.  Thus,  zero  tangential  electric  field 
boundary  conditions  near  this  third  interface  do  not  allow  the  electric  field  distribution  to  be 
“turned”  by  90°  that  was  detected  in  the  vertical  configuration  (Fig.  8). 

Concluding  it  is  worth  noting  that  the  position  of  the  resonance  in  the  frequency  scale 
for  the  horizontal  configuration  can  be  easily  shifted  by  changing  either  the  length  or  the 
height  of  the  obstacle  while  the  average  electric  field  in  it  can  be  adjusted  by  changing 
specific  resistance  of  the  SE. 

6.2  Frequency  response  optimization 

During  our  investigation,  two  possible  scenarios  for  optimizing  the  sensor  were 
considered.  First,  short  sensors  were  investigated  since  they  shouldn’t  introduce  too  much  of 
electric  field  distortion  in  the  waveguide.  By  changing  the  specific  resistance  of  the  SE,  it 
should  be  possible  to  exploit  the  interplay  between  conduction  and  displacement  currents  and 
obtain  the  desired  electric  field  ratio  for  different  frequencies.  Results  of  these  simulations  in 
cylindrical  coordinate  system  are  summarized  in  Fig.  20. 

As  one  can  clearly  see  from  the  figure,  it  is  possible  to  obtain  electric  field  ratios  close 
to  the  desired  1.23,  however  several  practical  problems  appear.  In  the  case  of  20  £Tcm 
material  (b)  it  is  possible  to  simultaneously  have  lower  than  50  Q  resistance  and  the  required 
ratio.  However,  small  change  of  the  dimensions  lead  to  the  substantial  increase  of  this  ratio, 
for  example,  the  change  of  d  from  0.4  to  0.8  mm  leads  to  an  increase  of  the  electric  field  ratio 
from  1.29  to  1.54.  Therefore,  the  practical  realization  of  such  a  sensor  could  face  some 
difficulties  due  to  unintentional  size  and  positioning  variations. 

In  addition  to  the  investigation  of  short  sensors,  we  also  performed  the  series  of 
calculations  for  longer  ones  at  fixed  d  =  1.77  mm  and  the  same  fixed  height  of  h  =  2  mm. 
These  results  for  both  considered  specific  resistances  are  also  presented  in  Fig.  20. 
Unfortunately,  only  the  monotonous  decrease  of  the  electric  field  ratio  was  obtained  for 
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Fig.  20.  The  ratio  of  the  average  electric  field  (E^)  calculated  at/=  14.5  GHz  to  (E^)  calculated  at/=  10.3  GHz 
versus  longitudinal  and  transverse  dimensions  of  the  SE  for  p=  10  Qcm  (a)  and  p=  20  Qcm  (b) 
material  and  fixed  h  =  2  mm  height.  The  number  near  the  point  corresponds  to  the  value  of  the  ratio  for 
the  particular  pair  of  d  and  /.  The  solid  curve  in  the  figure  corresponds  to  d  and  l  at  which  the  DC 
resistance  of  the  RS  is  50  Q.  Therefore,  above  this  curve,  the  resistance  of  the  RS  is  more  than  50  Q; 
below  this  curve,  it  is  less  than  50  Q.  Waveguide  radius  a  =10  mm. 
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Fig.  21.  The  ratio  (Eg,)(14.5)/(Eg,)(10.3)  versus  longitudinal 
and  transverse  dimensions  of  the  SE  for  p=  20  Qcm 
material  and  fixed  height  h  =  4  mm.  The  number  near 
the  point  corresponds  to  the  value  of  the  ratio  for  the 
particular  pair  of  d  and  /.  The  solid  curve  in  the  figure 
corresponds  to  d  and  /  at  which  the  DC  resistance  of 
the  RS  is  50  Q,  a  =  10  mm. 


considered  specific  resistances.  Such 
results  indicate  that  the  larger  sensors 
should  be  investigated  due  to  two 
reasons:  (i)  once  the  length  and/or 
height  of  the  SE  is  increased  the  lower 
resistance  is  obtained,  (ii)  the  larger 
size  of  the  SE  enhances  possibilities  to 
exploit  electric  field  increase  near 
resonant  frequencies. 

Results  of  the  simulations  for 
“taller”  SE  are  presented  in  Fig.  21. 
From  the  line  of  ratios  corresponding 
to  d-  1.57  mm  one  can  clearly  see, 
that  the  increased  height  leads  to  the 
increase  of  average  electric  field  ratio 
once  the  length  of  the  SE  reaches  3  - 
4  mm.  This  is  an  indication  that  higher 
sensors  become  “massive”  enough  to 
have  resonant  frequencies  near  the  end 
of  the  investigated  range  of  frequency 
and  therefore  more  suitable 


dependence  of  the  electric  field  on  frequency. 

Another  favorable  feature  of  this  layout  becomes  noticeable  when  exploring  ratios  for 
fixed  /  and  various  d  values  -  electric  field  at  higher  frequencies  tends  to  grow  faster  once  the 
width  d  is  reduced.  This  feature  is  particularly  important,  since  reduced  width  can  potentially 
reduce  the  reflection  coefficient  from  the  sensor.  It  also  should  be  emphasized,  that  the 
electric  field  ratio  changes  smoothly  once  dimensions  are  changing  within  0.5  mm  and 
therefore  such  a  sensor  should  be  more  suitable  for  the  reliable  production. 


6.3  Properties  of  the  optimal  sensor 

A  sensor  with  dimensions  hxdxl  =  4x0.4x3.5  mm3,  p=  20  Q-cm,  R  =  5.6  Q  was 
selected  for  a  more  detailed  investigation.  The  calculations  were  performed  at  more 


Fig.  22.  Dependences  of  the  average  electric  field  within  the  sensor  (a)  and  VSWR  (b)  on  the  angle  between 
center  of  the  SE  and  polarization  plane  of  the  incident  wave  for  several  frequencies.  The  ratio  Eg, I  sin  cp 
is  also  plotted  in  (a)  for/=  14.5  GHz. 
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b)  hxdxi  =  4x0,4*3,5  mm3,  p  =  20  Q  cm 


f,  GHz 


Fig.  23.  Dependences  of  the  sensitivity  (a)  and  normalized  sensitivity  (b)  on  frequency  for  the  optimal  sensor  in 
the  horizontal  configuration  for  several  different  polarization  angles  indicated  in  the  figure. 


frequencies  and  dependencies  on  the  polarization  angle  were  considered. 

Dependencies  of  the  averaged  electric  field  and  voltage  standing  wave  ratio  on 
polarization  angle  for  three  frequencies  are  presented  in  Fig.  22.  From  Fig.  22a  follows  that 
the  angular  dependence  of  the  average  electric  field  in  the  SE  at  14.5  GHz  almost  perfectly 
matches  the  sin  cp  dependence  that  is  characteristic  of  the  regular  E ^  component  of  the 
incident  Hu  wave.  It  seems  that  similar  angular  dependence  is  typical  of  the  other  frequencies 
as  well.  Such  property  can  be  attributed  to  the  geometry  of  the  SE.  It  is  thin  in  (p  direction  and 
its  metal  contacts  are  perpendicular  to  E '9  electric  field  component.  It  is  seen  that  the  SE 
perfectly  detects  the  maximum  of  the  E '9  amplitude  {cp-  90°)  at  every  frequency  therefore  it 
means  that  such  device  can  act  as  a  nearly  perfect  sensor  for  the  measurement  of  a  local  E 9 
value.  Unfortunately,  metal  contacts  perpendicular  to  the  walls  of  the  waveguide  introduce 
substantial  reflections  once  the  SE  is  placed  in  the  position  where  E,  value  is  at  its  maximum. 
As  follows  from  Fig.  22b,  the  VSWR  increases  substantially  for  small  preaching  almost  the 
value  of  2.2  when  the  SE  is  placed  in 
the  position  cp  -  0. 

As  in  a  case  of  the  vertical 
configuration,  the  current  sensor  also 
should  be  rotatable,  but  the  maximum 
of  the  signal  should  be  found.  In  such  a 
situation,  the  SE  is  perpendicular  to  the 
polarization  plane  of  the  incident  wave.  £: 

As  follows  from  Fig.  22b,  the  % 
reflection  from  the  SE  can  also  indicate 
the  desirable  position  of  the  SE  in 
respect  of  polarization  plane.  At 
cp  -  90°  VSWR  is  at  its  minimum. 

An  extensive  investigation  has 
been  performed  to  reveal  more  details 


hxdxl  -  4x04x3.5  mm3,  p  =  20  Q  cm 


of  the  dependence  of  the  sensitivity  on 
frequency.  Results  of  the  simulation 
are  presented  in  Fig.  23.  As  it  can  be 
expected,  sensitivity  increases  with  the 
increase  of  the  angle  cp.  It  was  already 


Calculated  dependence  of  the  relative  resistance 
change  of  the  optimal  RS  on  power  transmitted 
through  the  waveguide  at  f-  13.6  GH.  Solid  line 
shows  results  calculated  using  (16),  dotted  line 
corresponds  to  linear  dependence. 
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mentioned  that  the  increase  of  the  angle  leads  to  the  growth  of  the  component  E ^  of  the 
incident  wave  near  the  SE.  Unfortunately,  the  absolute  value  of  the  sensitivity  varies  within 
±14  %  over  the  entire  frequency  range.  However,  normalized  sensitivity  remains  the  same 
within  1.2  %  for  every  fixed  frequency,  confirming,  that  this  device  can  be  of  some  use  as 
nearly  the  ideal  sensor  for  a  fixed  frequency  to  determine  the  position  of  the  polarization 
plane  in  the  waveguide. 

Once  the  electric  field  strength  in  the  SE  is  determined,  the  dependence  of  the  relative 
resistance  change  on  a  power  transmitted  through  the  waveguide  can  be  calculated  in  the 
similar  manner  as  in  subsection  5.2,  Fig.  lib.  Calculation  results  are  shown  in  Fig.  24  by  a 
solid  line.  Dotted  line  corresponds  to  the  linear  dependence  of  A R/R  on  P  that  is  characteristic 
of  the  warm  electron  region.  It  is  seen  that  at  a  maximum  power  roughly  35%  relative 
resistance  change  of  the  SE  is  observed. 

Due  to  a  small  width  of  the  optimal  SE,  the  averaged  electric  field  in  it  consists  of  E ^ 
component  only.  The  amplitude  of  the  component  Er  is  roughly  10"4  and  Ez  is  even  smaller. 
The  distributions  of  the  electric  field  components  in  the  waveguide  section  with  the  SE  are 
shown  in  Fig.  25.  The  distribution  of  the  components  E ^  (a)  and  (d),  Er  (b)  and  Ez  (c)  in 
different  (pOz  planes  are  shown.  It  is  seen  that  due  to  a  small  width  of  the  SE,  the  perturbation 
of  the  regular  components  of  Hu  wave  at  cp  =  90°  is  small  (a  and  b).  VSWR  for  this  particular 
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Fig.  25.  Distributions  of  the  components  of  electric  field  in  the  waveguide  with  the  SE  on  the  different  surfaces 
q$) z  at/=  10.3  GHz:  at  cp-  90°  (a)  and  cp  =  0°  (d)  for  r  =  (/-l)Ar  ;  Er  at  cp=  90°  (b)  for  r  =  (/-0.5)Ar; 
Ez  at  #?=90°  (c)  for  r=(z-0.5)Ar;  Calculation  parameters:  Ar  =0.025,  Az  =0.05,  A<^=2.8°, 
At  =  6.10-10'4.  Dimensions  of  the  SE  are  hxdxl  =  4.0x0.4x3.55  mm3,  p=  20  Q-cm,  a  =  10  mm. 
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case  is  roughly  1.02.  Some  perturbation  of  the  regular  wave  is  seen  in  the  planes  close  to  the 
SE.  The  component  Ez  that  is  absent  in  the  regular  wave  is  excited  but  its  amplitude  is  very 
small  less  than  1%  in  the  center  of  the  waveguide.  It  becomes  more  significant  in  the  obstacle 
or  near  it.  The  distribution  of  electric  field  changes  significantly  when  the  position  of  the  SE 
coincides  with  polarization  plane  (d).  Partly  standing  wave  between  the  SE  and  the  excitation 
plane  appears  demonstrating  strong  reflection  from  the  SE.  VSWR  for  this  particular  case  is 
roughly  1.6. 


6.4  Simulations  in  Cartesian  coordinate  system 


Semiconductor  technology  usually  deals  with 
parallelepiped  shaped  objects  cut  from  plates  with 
parallel  surfaces  (Fig.  26a).  The  SE  considered  in  the 
cylindrical  coordinate  system  has  tapered  cross  sectional 
shape  (Fig.  26b).  Variation  of  thickness  becomes  more 
substantial  when  the  height  of  the  sensor  approaches  the 
radius  of  the  waveguide.  Wishing  to  elucidate  the 
influence  of  the  shape  to  the  sensors  sensitivity,  we 
performed  several  FDTD  calculations  using  Cartesian 
grid  and  the  parallelepiped  shaped  SE  as  it  is  shown  in 
Fig.  26a. 

A  comparison  of  the  sensitivity  of  the  SE  on 
frequency  for  two  coordinate  systems  is  presented  in  Fig. 
27.  As  one  can  clearly  see,  sensitivities  coincide  in  a  high 
frequency  region  demonstrating  up  to  1 1  %  difference  in 
a  low  frequency  region.  The  obtained  discrepancy  might 
be  attributed  to  the  difference  in  shapes. 


Fig.  26.  Sensor’s  shapes  fitting  Carte¬ 
sian  (a)  and  cylindrical  (b) 
coordinates. 


hxdxl  =  4x0. 5x3.5  mm  ,  /?-  20  Q  cm 


6.5  Conclusions 

Two  possible  scenarios  for 
obtaining  a  nearly  flat  frequency 
response  of  the  horizontal 
configuration  sensor  (contacts  on  the 
sidewalls)  were  considered.  One 
possibility  is  to  build  very  short, 
narrow  and  relatively  high  sensor. 

However,  in  this  case,  its  resistance  is 
over  50  Q  and  its  sensitivity  is  very 
susceptible  to  small  changes  of 
dimensions.  Another  possibility  is  a 
narrow  plate  shaped  sensor.  The 
optimal  SE  found  at  a  moment  should 
have  the  following  parameters: 
hxdxl  =  4x0.4x3.5  mm3,  p=  20  Qcm, 

R  =  5.6  Q.  Such  sensor  retains  the 
dependence  of  its  normalized 

sensitivity  on  frequency  for  changing  the  angle  between  sensor’s  position  and  polarization 
plane  of  the  incident  wave,  but  the  sensitivity  changes  quite  substantially  ±14  %  of  the 
average  value  in  the  frequency  range  10.3  -  14.5  GHz.  Therefore,  such  sensor  would  be 
applicable  only  for  fixed  frequency  but  variable  polarization  experiments.  In  general  vertical 


Fig.  27.  Dependences  of  the  sensitivity  on  frequency  for  the 
sensor  computed  in  two  coordinate  systems  indicated 
in  the  figure. 
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configuration  RS  and  especially  the  array  of  two  sensors  presented  in  subsection  5.3  is  more 

preferable. 

7  Conclusions 

1 .  Making  use  of  a  FDTD  method,  the  peculiarities  of  the  interaction  of  the 
semiconductor  obstacle  placed  on  a  wall  of  the  circular  waveguide  with  a  mode  Hu 
(TEn)  have  been  investigated.  Two  different  structures  have  been  considered  with  the 
contacts  on  top  and  bottom  planes  (vertical  configuration)  and  with  the  contacts  on 
sidewalls  (horizontal  configuration).  Both  these  configurations  can  be  considered  as  a 
prototypes  of  the  RS  devoted  for  HPM  pulse  measurement  in  a  circular  waveguide. 

2.  A  program  based  on  the  FDTD  method  was  written,  tested,  and  used  for  the 
calculation  of  the  electromagnetic  field  components  in  the  circular  waveguide  with  the 
semiconductor  obstacle  exposed  by  Hu  mode  regular  wave. 

3.  Considering  the  sensitivity  of  the  RS  in  the  circular  waveguide,  it  was  shown  that  the 
decrease  of  the  sensitivity  with  frequency  occurs  due  to  the  wave  dispersion  in  the 
waveguide.  To  compensate  it  the  increase  of  the  average  electric  field  in  the  sensing 
element  of  the  RS  with  frequency  is  required. 

4.  Considering  insulating  obstacle  in  the  circular  waveguide,  it  was  found  that  the 
resonances  occur  in  the  structure  when  the  length  of  the  obstacle  becomes  equal  to  the 
integer  number  of  Ad/2,  where  Ad  is  the  wavelength  of  microwaves  in  the  structure.  In 
general,  the  resonance  type  and  frequency  depends  on  the  position  of  the  contacts 
therefore  the  resonances  are  different  for  the  vertical  and  horizontal  configurations. 

5.  The  resonance  in  both  configurations  can  be  shifted  to  the  lower  frequency  by 
increasing  the  length  of  the  SE;  the  averaged  electric  field  amplitude  in  it  decreases 
when  the  specific  conductivity  of  the  SE  grows.  For  the  vertical  configuration,  the 
increase  of  the  width  of  the  SE  acts  on  the  resonance  position  in  the  same  way  as  the 
length  while  the  change  of  the  height  of  the  SE  has  little  influence  on  the  resonance 
frequency.  For  the  horizontal  configuration,  on  the  contrary,  the  resonance  frequency 
can  be  shifted  to  the  lower  frequency  by  increasing  its  height  while  the  change  on  the 
width  has  little  influence  on  the  resonance  position. 

6.  The  optimal  vertical  configuration  RS  found  at  a  moment  should  have  the  following 
parameters:  p/a  =  20  Q,  hlaxdlaxlla  =  0.05x0.38x0.5.  Calculated  sensitivity  variation 
within  frequency  range  fjf =0.85-0.60  for  the  optimal  SE  placed  in  the  polarization 
plane  of  the  incident  wave  is  roughly  ±3.5%,  where  fc  is  the  cutoff  frequency  for  Hu 
mode.  The  DC  resistance  of  the  RS  R  =  21  Q,  calculated  VSWR  slightly  increases 
with  frequency,  but  it  is  less  than  1.03  within  considered  frequency  range.  Since  the 
output  signal  from  the  RS  depends  on  its  position  in  respect  to  the  polarization  plane, 
the  rotatable  SE  should  be  used  to  find  the  minimum  signal  from  the  RS.  In  this  case, 
polarization  plane  coincides  with  the  position  of  the  center  of  the  SE. 

7.  The  alternative  RS  has  been  proposed  consisting  of  two  SE  with  parameters: 
p/a  =  20  Q,  hlaxdlaxlla  =  0.05x0.29x0.5  placed  at  90°  from  each  other.  It  was  shown 
that  the  sum  of  their  signals  within  ±7%  is  independent  of  the  position  of  the 
polarization  plane  in  the  waveguide.  The  DC  resistance  of  each  RS  R  =  24  Q, 
calculated  value  of  the  VSWR  for  two  SE  was  less  than  1.05  within  considered 
frequency  range. 

8.  The  optimal  horizontal  configuration  RS  found  at  a  moment  should  have  the  following 
parameters:  p/a  =  20  Q,  hlaxdlaxlla  =  0.4x0.04x0.35.  Calculated  sensitivity  variation 
within  frequency  range  fjf  =0.85-0.6  for  the  optimal  SE  placed  in  the  plane 
perpendicular  to  the  polarization  plane  is  roughly  ±14%.  The  DC  resistance  of  the  RS 
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R-  5.6  CX  Since  the  output  signal  from  the  RS  depends  on  its  position  in  respect  to  the 
polarization  plane,  the  rotatable  SE  should  be  used  to  find  the  maximum  signal  from 
the  RS.  In  this  case,  polarization  plane  is  perpendicular  to  the  position  of  the  center  of 
the  SE.  Calculated  VSWR  strongly  depends  on  polarization  plane  position,  gaining 
minimum  value  at  a  maximum  signal. 

9.  Comparing  optimal  sensors  found  for  both  configurations,  it  is  seen  that  both  could 
measure  HPM  pulse  power  up  to  the  breakdown  in  the  waveguide.  It  seems  that 
vertical  configuration  RS  and  especially  the  array  of  two  sensors  is  more  preferable. 
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1  Appendix 

Written  programs  have  been  tested  by  comparing  computed  results  with  the  analytical 
solution  or  results  obtained  using  the  other  programs.  In  the  subsection,  the  tests  are 
described. 

1.1  Comparison  with  analytical  solution 

It  worthwhile  to  mention  that  even  small  inaccuracies  those  are  left  in  the  FDTD 
program  leads  to  unstable  solution:  small  inaccuracies  are  multiplied  in  the  successive  steps 
when  the  new  values  of  the  components  of  the  electromagnetic  field  are  updated  from  the  old 
ones  leading  to  the  infinite  growth  of  them.  Therefore,  the  preliminary  testing  of  the  program 
is  very  simple.  One  should  get  non-infinite  values  of  the  electromagnetic  field  components 
after  a  few  periods  of  excitation  of  the  regular  wave  in  the  waveguide.  The  next  step  that  we 
undertook  in  program  testing  was  the  calculation  of  the  reflection  coefficient  from  the  width  d 
dielectric  wafer  that  is  tightly  inserted  into  the  circular  waveguide  so  that  it  fully  covers  its 
window.  From  the  one  hand,  the  reflection  coefficient  from  such  obstacle  can  be  calculated 
using  the  program  described  in  Section  4.  From  the  other  hand,  it  can  be  computed 
analytically  considering  plane  wave  incident  at  angle  6  from  the  normal  to  the  dielectric  plate 
of  width  d.  For  Hu  mode  the  wave  with  E  in  plane  of  incidence  should  be  considered. 

a)  b) 


Fig.  28.  Dependences  of  the  reflection  (a)  and  transmission  coefficients  (power)  on  the  width  of  the  dielectric 
wafer.  The  wafer  is  inserted  into  circular  waveguide,  mode  Hu.  Lines  show  analytical  solution,  points 
demonstrate  calculations  results  of  the  FDTD  program. 

Calculation  results  obtained  using  our  FDTD  program  and  the  analytical  solution  for 
two  different  frequencies  are  shown  in  Fig.  28.  The  radius  of  the  waveguide,  the  relative 
dielectric  constant  and  the  electromagnetic  wave  frequency  are  given  in  the  figure.  It  is  seen 
that  at  some  thickness  of  the  wafer  the  Fabri-Perot  resonance  conditions  are  fulfilled,  the 
wafer  brightens  and  the  reflection  coefficient  goes  to  zero.  Good  coincidence  between  the 
numerical  results  and  analytical  solution  supports  the  validity  of  the  FDTD  program  used  in 
the  calculations. 


